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ABSTRACT 


The finned tube contactor of a liquid desiccant (triethylene 
glycol) dehumidifier system is modeled with a fifth order differential 
equation to simulate the simultaneous heat and mass transfer associated 
with conditioning the air. Potential energy saving aspects of the sys- 
tem are discussed. Computer simulation results presented on psychro- 
metric charts show that significant performance is possible with contac- 
tors of one half the recommended size. The cycle may be regenerated at 
120 to 140°F depending on the system design requirements. 


Figure 4 may be used to estimate system performance. 
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NOMENCLATURE AND UNITS 


Abbreviations 
LDD Liquid Desiccant Dehumidification 
TEG Triethylene Glycol 


English Letter S ols 


ao Cooling area per unit volume for heat transfer, 
fever 

ay Interfacial area per volume for heat transfer, 
Bee 

ay Interfacial area per volume for mass transfer, 
fest? 

Go€ Heat capacity at constant pressure, BTU/1b 

Cy Humid heat, BTU (of water vapor-gas mixture) / 
(1b dry gas)(°F) defined by Eq. (A.22) 

Co Humid heat, BTU (of TEG-water solution) /(1b TEG) 
(°F) defined by Eq. (A. 33) 

Oe Root diameter of tube, ft 

dD Diffusion coefficient for water vapor in air, 
ft /he 

Fo Mass Transfer coefficient, 1b water/(hr) (es) 
defined by Eq. (C.1) 

g Acceleration due to gravity, elie 

G Superficial gas (dry air plus water vapor) mass 


velocity, yin yCee.) 





max 


Gas mass velocity at minimum cross section, 

1b/ (hr) (£t”) 

Superficial solvent (dry air) mass velocity, 

1b/ (hr) (£t”) 

Contact coefficient of heat transfer between gas 
and liquid, BTU/ (hr) (£t2) (°F) defined by Eq. (C.2) 
Contact coefficient of heat transfer between gas 
and liquid corrected for simultaneous mass trans- 
fer BTU/ (hr) (£t-) (°F) defined by Eq. (A.15) and 
(A.16) 

Contact coefficient of heat transfer between 
liquid and cooling fin, BTU/ (hr) (ft) (°F) defined 
by Eq. (C.13) 

Conduction coefficient of heat transfer in liquid, 
BTU/ (hr) (ft) (°F) defined by Eq. (C.11) 

Enthalpy, BTU/1b 

Heat of dilution (mixing), BTU/1b 

Thermal conductivity, BTU/ (hr) (ft) (°F) 

Effective fin length, in. 

Superficial liquid (TEG plus water) mass velocity, 
1b/ (hr) (£t”) 

Superficial solvent (TEG) mass velocity, 


ih ne) (ec) 






pypreey rviiii 


sh 


SV 


Mass transfer flux of water vapor, iby Chey (Et) 
defined by Eq. (A.14) 

Partial pressure in equilibrium with a liquid 
solution, mm Hg 

Total pressure, mm Hg 

Fin pitch, fins/in; Chapter II: pressure 

Heat removed by cooling coil, BTU/hr defined by 
Eq. (A.28) 

Exothermic heat of dilution, BTU/hr defined by 
Eq. (A.30) 

Heat transfer flux for sensible heat, BTU/ (hr) 
(£t7) 

Chapter II: universal gas constant; 

Appendix B: bond graph resistance 

Spacing between fins, in. 

Interfacial surface, aoe cross section 
Thickness of aluminum cooling fin, in; tempera- 
ture, °F 

Chapter II, temperature, °F 

Horizontal tube spacing, in. 

Vertical tube spacing, in. 

Over-all coefficient of heat transfer between the 
liquid and cooling water, BTU/(hr) (ft) (°F) 


Water velocity through tube, ft/sec 






Pee TPT eeteaty 


y* 


Greek Letter Symbols 
It 


5 


10 


Mass fraction concentration in the bulk liquid, 
1b water/(1b water + 1b TEG) 

Mass ratio concentration in the bulk liquid, 

1b water/1b TEG 

Mass fraction concentration in the bulk gas, 

1b water vapor/(1b water vapor + lb dry air) 
Mole fraction concentration in the gas in equilib- 
rium with the bulk liquid, moles water vapor/ 
(moles water vapor + moles dry air) 

Mole fraction concentration in the bulk gas, 
moles water vapor/(moles water vapor + moles dry 
air) 

Mass ratio concentration in the bulk gas, 1b 
water vapor/(1b dry air) 

Height of finned tube contactor (parallel to 


flow direction), ft. 


Liquid loading per unit width of fin, 1lb/(hr) (ft) 
TEG solution thickness, ft. 

Fin efficiency, dimensionless 

Latent heat of vaporization, BTU/1b 

Viscosity at bulk temperature, 1b/(hr) (ft) 
Density, yee 


Ratio of free flow area to face area, dimension- 


less 





v 


Dimensionless Groups 


Pr 


se 


Subscripts 
A 


AL 


£56 


ab 


Bond Graph Symbols 


ial 


Chapter II: specific volume, £t°/1b 


Prandtl number (ic 7k) 


Schmidt number (/p dD) 


Substance A, water vapor 

Substance A, liquid water 

Substance B, dry air 

Substance C, TEG 

Substance D, water plus TEG solution; dilution 
Cooling water 

Gas 

Interface 

Liquid 


At the reference temperature 


Capacitance 

Mass capacitance 
Thermal capacitance 
Mass transfer flux 
Heat flux 

Resistance 
Diffusion resistance 


Interfacial resistance to mass transfer 
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Reaction resistance 

Heat of dilution flow source 

Latent heat flow source 

Temperature 

Transformer 

Parallel junction, flows sum to zero 
Series junction, efforts sum to zero 


Power bond 


Signal bond 
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I. INTRODUCTION 


Definitions 
Liquid desiccant dehumidification is the reduction of the water 

vapor content of a given volume of air by a liquid sorbent which has the 
property of extracting and holding water vapor brought into contact with 
it. The terms sorbent and desiccant refer to materials with a large ca- 
pacity for moisture compared to their volume and weight. The amount of 
water held by a desiccant will increase or decrease depending on whether 
the vapor pressure of water above the desiccant is less than or greater 


than the partial pressure of water vapor in the air in contact with it. 


Applications 


Dehumidification, whether by liquid desiccants or by other means, 
has important commercial applications. These include preserving ships 
and other surplus equipment which would otherwise deteriorate, and pre- 
venting condensation damage to cargo and internal structures of ships' 
holds which transport dry cargo. 

Liquid desiccant dehumidification (LDD) was investigated in the 
1950's and found to be operational [7],[22]. Currently LDD equipment 
designed for low pressure steam regeneration is commercially available 
and is used in the drying of natural gas prior to pipeline transporta- 
tion. Both hospitals and the food processing industry use LDD not only 
for temperature and humidity control but also for the bactericidal prop- 
erties of the family of glycol desiccants. Accepted high equipment cost 


prevents wider application of LDD for the general purpose of space de- 
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humidification or space cooling. These costs are on the order of $10 
per square foot of contactor in typical LDD systems which remove about 
500 pounds of water per hour with dehumidification and regeneration con- 


tactor areas of 5000 sq. ft. each [27]. 


Potential Savings 


If the potential to regenerate an LDD cycle at low temperatures is 
exploited and waste heat or solar energy is utilized, then potential 
energy savings can be applied against acquisition cost. Such cost may 
be minimized by rational rather than empirical equipment selection. 

A review of current space dehumidification and cooling practice 
will indicate the potential energy savings of LDD. 

In conventional air conditioning systems the air is cooled signifi- 
cantly below the desired temperature to effect moisture removal. (Refer 


to Figure 1). 


CHILLED WATER ENERGY 


COOLING | VERY COOL, DRY AIR 
COIL 


Figure 1. Conventional Air Conditioner with Reheat 






REHEAT 
COIL 





WARM, MOIST AIR 








15 


Additional energy is used as this very cool, dry air is then reheated 
to the desired space temperature. The popularity of this system stems 
from the previously low cost of energy and the simple control system re- 
quired. One dependent variable, outlet humidity ratio, is uniquely con- 
trolled by the chilled water temperature through the condensation of ex- 
cess moisture above the dew point temperature. The other dependent vari- 
able, desired space temperature, is now uniquely controlled at a constant 
humidity by the quantity of reheat energy added. 

The control system complexity increases when an attempt is made to 
Save reheat energy by mixing the very cool dry air with warm moist air 


bypassed around the cooling coil. (Refer to Figure 2.) 


CHILLED WATER 


WARM, MOIST AIR | cOOLING| VERY COOL DRY AIR 
COIL 
BYPASS AIR 


Figure 2. Conventional Air Conditioner with Bypass Mixing 












COOL, LESS DRY AIR 





Both of the two independent variables, chilled water temperature and bypass 
air flow rate, simultaneously modify both of the dependent variables, 
outlet air temperature and humidity ratio. 

The LDD system can save the energy of excess chilling for moisture 
removal through its ability to dehumidify at ambient temperatures. Re- 
moval of the latent heat of condensation is the primary function of the 


cooling water. (Refer to Figure 3.) 
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CONCENTRATED 


LIQUID DESICCANT COOLING WATER 


AIR—LIQUID | COOL DRY AIR 
CONTACTOR 


WARM MOIST AIR 
DILUTE 


LIQUID DESICCANT 






Figure 3. Liquid Desiccant Dehumidification System 


The dehumidification contactor of the LDD system, a three fluid heat ex- 
changer in which mass transfer between two of the fluids occurs, can 
supply cool, dry air without further heating or cooling. Thus, the re- 
heat energy can also be saved. 

In the LDD contactor the two dependent variables, air temperature 
and humidity ratio, are controlled by the independent variables, liquid 
desiccant concentration and cooling water temperature. A tradeoff for 
the potential energy savings of excess chilling and unnecessary reheat 
is a control system of sophistication similar to that of the bypass sys- 
tem shown in Figure 2. 

The LDD cycle compares favorably to an alternative technique for 
space cooling such as the lithium bromide absorption cycle. While the 
absorption cycle is closed and the desiccant cycle is open, both systems 
have coefficients of performance on the order of 0.6 [11]. The absorp- 


tion cycle requires regeneration temperatures around 200°F, while the 
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LDD cycle can be regenerated at temperatures around 140°F. Of these 

two systems the absorption system has been the subject of considerable 
analysis and simulations [9]; however, there is no well-established tech- 
nique for space cooling or dehumidification with low temperature regener- 


ation. 


Purpose 


In view of the above, a model for analysis of the liquid desiccant 
dehumidifier contactor of Figure 3 is presented as the basis for a ra- 
tional design. The primary function of the model is to assure that the 
contactors are not over-specified so that acquisition costs may be kept 
as low as possible. The contactors (dehumidification and regeneration) 
represent a significant portion of the total cost of an LDD system. 

Empirical equipment selection of the contactors is standard prac- 
tice and can result in overspecification to ensure proper operation. 
Manufacturers’ application recommendations specify equipment which oper- 
ates with the exit conditions of air and liquid temperatures, humidity 
ratio, and desiccant concentration at equilibrium. 

By contrast, rational design implies more than black box modeling 
of the contractor. The computer model of the contactor which was devel- 
oped is based on a complete set of energy balances, material balances, 
rate equations, and equilibrium relationships. Through computer simu- 
lation of the heat and mass transfer in the contactor, the rate at which 
the temperatures, humidity ratio, and desiccant concentration approach 
equilibrium as a function of flow distance can be studied. With a speci- 


fic design the contactor flow length and associated cost can be compared 
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to the incremental heat and mass transfer of additional contactor flow 


length increments. 





wg) 


IIT. PSYCHROMETRY 


The Psychrometric Chart 


Here psychrometry is concerned with the determination of the prop- 
erties of the mixture of air and water vapor as the mixture flows through 
the contactor. These properties are conveniently diagrammed in an abbre- 
viated psychrometric chart, Figure 4, which plots dry bulb temperature, 
humidity ratio, wet bulb temperature, saturation temperature, and rela- 
tive humidity for the air-water vapor mixture. 

Dry bulb temperature represents the air and liquid contact tempera- 
ture and is the abscissa. 

Humidity ratio (H.R.) is the pounds of water carried by one pound 
of dry air. Alternatively, grains per pound of dry air (grains) is used 
where 7000 grains equal one pound of water. The ordinate shows both 
humidity ratio and grains. 

Relative humidity is the ratio of partial pressure of the water 
vapor to the vapor pressure of water at the air temperature and is indi- 
cated by the dotted lines. Relative humidity is expressed in percent. 

Wet bulb temperature lines run upward and to the left to the satu- 
ration curve and very nearly represent the path followed by adiabatic 
(spray) cooling. 

Dew point temperature is the temperature to which the mixture must 
be cooled to become saturated and is found by following a horizontal line 
to the saturation temperature curve. 


Saturation temperature represents the maximum moisture the air can 
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DRY BULB TEMPERATURE, Deg. F 


Triethylene Glycol Equilibrium Humectant Values on Psychrometric Chart 


Figure 4, 
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hold before fogging ensues. 

Enthalpy and specific volumes are plotted on more detailed charts 
[26]. Any two of these parameters specify the air-water vapor mixture 
and allow all the others to be determined. The remaining parameter, 
humectant value, is discussed later under its own section. 

Dry bulb temperature and humidity ratio were chosen as the two de- 
pendent variables of interest because: (1) dry bulb temperature is 
easily measured and has significance to people; and (2) the humidity 
ratio directly represents the water content of the air and is based on 


one pound of dry air which is a constant in any psychrometric process. 


Psychrometric Chart Processes 


The process of heating or cooling moist air above the dew point is 
represented by a horizontal line since the humidity ratio is unchanged. 
This process is known as sensible heating or cooling. 

The process of cooling below the dew point entails condensation of 
vapor and is represented by a horizontal line extending left to the satu- 
ration curve where the process follows the curve to a lower humidity 
ratio. This is roughly the dehumidification process of conventional air 
conditioning. 

A downward vertical line on the chart represents an isothermal dehu- 
midification process or latent cooling. This may be accomplished directly 
by the LDD system with cooling to remove only the latent heat of conden- 


sation. 


Liquid Desiccants 


A liquid desiccant removes water vapor from the air at a rate direc- 
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tly proportional to the difference between the partial pressure of the 
water vapor immediately above the liquid desiccant. The partial pres- 
sure of water vapor in the air is reduced and the liquid desiccant is 
diluted in the process until equilibrium occurs. The process may be re- 
versed by heating the desiccant. That is, the liquid desiccant is con- 
centrated as the partial pressure of water vapor in the air increases. 
The vapor pressure of the liquid desiccant is a function of temperature 


and concentration. 


Humectant Values 

Equilibrium humectant values are the percent by weight of a liquid 
desiccant in water solution required to maintain vapor pressure equili- 
brium in contact with air of various temperatures and humidity ratios. 
The humectant values for triethylene glycol (TEG) were added to the 
standard psychrometric chart to produce a convenient diagram revealing 
the relationship between TEG concentration and the air-water vapor mix- 
ture at equilibrium. The humectant values on Figure 4 were developed 
from graphs of the vapor pressure of water over triethylene glycol of 
various weight concentrations and temperatures [28] [37]. These equilib- 
rium humectant values give an indication of the TEG concentration necess- 
ary to dehumidify to a specified outlet condition with a constant TEG 
concentration. 

The solid lines which are the TEG equilibrium humectant values rep- 
resent slightly different TEG weight concentrations in the 90% to 99% 
range depending on the data base: Union Carbide [37] or Dow Chemical [281] 


A third source [12], presented data for only three concentrations, 40%, 
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60%, and 80% and these were consistent with the first two sources within 
9%. The vapor pressures presented by Union Carbide and Dow Chemical 


differ by approximately 20% at 90% concentration and 100% at 98% concen- 
tration. Union Carbide data presented the higher vapor pressures. The 
more conservative humectant values during dehumidification are derived 
from the Union Carbide data while the more conservative humectant values 
during regeneration are derived from the Dow Chemical data, 

In the area of Figure 4 in which dehumidification takes place, 70 
to 80 °F and 0.010 H.R., the change in TEG humectant value with respect 
to temperature is about 1%/°F. At a 0.025 H.R. and 120 to 130 °F where 
regeneration can take place, the ratio is about 0.2%/°F which indicates 
that TEG regeneration should not be particularly sensitive to regener- 


ation temperatures. 


Vapor Pressure Prediction 


An equation for the equilibrium partial pressure above the liquid 
desiccant is necessary for the computer simulation. In addition to 
accuracy, a suitable equation will have continuous derivatives over its 
range and will be relatively efficient with respect to calculation time. 
Two relationships with these characteristics which predict vapor pres- 
sure as a function of TEG contact temperature and TEG weight concen- 
tration were developed with the aid of a curve fitting program and a 
programmable calculator: 


From Dow Chemical data, 


21.36(98-P_) 98-P 
p* [mm Hg] = exp ae + 20.5749 - ——-—* 5 + = + 
6259+ (98-P,,) 
p. 13.5 0.6024 
(2 = cs -1 (221) 
99 


0.28294 
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From Union Carbide data, 


—9259 1 1.5 


p* [mm Hg] = exp T + 460 + 20.407 - 1-9333590.5-P + 
/ Pr ,3)0.167 
1- (95-5) - 1 (232) 
Os26921 


where T = contact temperature, °F 


¥r 


TEG weight fraction concentration, %. 
Both relationships are accurate within about 3% within the ranges 
40 < T < 300, 75 < C < 99 and are slightly less accurate for 0 < C < ioe 


Eq. (2.2) was used exclusively in the simulation. 


Vapor Pressure Correlation 
The method used in synthesizing equations (2.1) and (2.2) is based 


on the Clausius—-Clapeyron equation which relates the slope of the vapor- 


pressure curve to the latent heat of vaporization 


gee r (2.3) 
dT T(V¢ - Vv) 

where YG = specific volume of the saturated vapor 

Vv, = specific volume of the saturated liquid 
A = latent heat of vaporization 

Using the ideal-gas law and neglecting Vy, in comparison with Ye 

dinp = = = Ad? (2.4) 
RT 


If X} is fairly constant over the temperature range 


_ A 
In P = = RT + constant (2.5) 
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Equation (2.5) suggests that a plot of log P against + will be linear 


over small temperature ranges. Antoine's equation published in 1888 


allows for slight variation of A with temperature [26] 


b 
it 4p ( 





InP =a- (2.6) 


and is probably the best vapor-pressure correlation equation for a simple, 
three-constant (a, b, and c) equation. Equation (2.6) is the basis for 


equations (2.1) and (2.2). 
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III. LIQUID DESICCANT DEHUMIDIFICATION SYSTEM 


Operating Cycle 


One system for liquid desiccant dehumidification (LDD) is shown 
in Figure 5. The dehumidifier section on the tere accepts warm moist 
air from the atmosphere and contacts the air with the liquid desiccant 
in the dehumidifier contactor. Intimate contact between the air and 
desiccant over a large area occurs through the formation of a thin film 
of desiccant as it is sprayed over the extended surface contactor. The 
heat of condensation and the heat of dilution flow through the liquid to 
the cooling water which can be sea water, well water, or water from a 
cooling tower. The quantity of cooling water required is strongly de- 
pendent on the sensible and latent heat removed from the air by the 
liquid desiccant and weakly dependent on the desiccant circulation rate 
between the dehumidifier and regenerator. 

The air leaves the contactor flowing downward with entrained desic- 
cant and loses most of the liquid to the sump while crossing it. The re- 
mainder of the entrained liquid is eliminated as the cool, dry air flows 
upward to the controlled space. 

The dry bulb temperature of this cool, dry air is weakly dependent 
on the water removed from the air but is strongly dependent on the temp- 
erature of the liquid desiccant. Accordingly, the air flowing through 
the dehumidifier can be warmed or cooled as it is dried depending on the 
temperature of the liquid desiccant and the temperature of the cooling 
water. 


Concentration of the liquid desiccant is more easily achieved with 
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dry air which would normally be exhausted from the dehumidified space, 
but can also be achieved with outside air. The heat of vaporization 

and dilution is supplied through the heating coil in the regenerator 
contactor to the liquid desiccant as it contacts the relatively dry air. 
Heating evaporates the liquid desiccant though at a much reduced rate 
than that of the water. A reflux coil maintained at about 90 °F is used 
as a condenser for the desiccant to prevent excessive loss in the hot 
exhaust. 

Figure 5 also shows the four pumps necessary for the liquid desic- 
cant movement--two spray pumps and two feed pumps. The liquid-to-liquid 
heat exchanger between the conditioner sump which is slightly warmer 
than the cooling water, and the concentrator sump which is slightly 
cooler than the heating water, improves the system thermal efficiency. 
An air-to-air heat exchanger between the hot wet exhaust and the cooler 
inlet air of the regenerator will also increase the system thermal effi- 
ciency. 

Figure 6 shows a similar system with a combined sump which requires 
only one pump. While this system is more compact and uses fewer pumps, 
thermal efficiency can be degraded since the sump operates at a mean 
temperature. With equal dehumidifier and regenerator desiccant recircu- 
lation rates the sump operates as a counter flow heat exchanger with an 
effectiveness of 0.5, or less than is normally expected from the counter 
flow heat exchanger of Figure 5. 

In either system the solution in the sump will require treatment as 


it will concentrate solids absorbed from the air. 
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Air-Desiccant Contactor 

The contactor is the part of the system in which the air is brought 
in contact with the liquid desiccant for dehumidification or regener- 
ation. The contactor should have high heat and mass transfer coeffici- 
ents but yet have relatively low air flow resistance. The contactor 
should also be light, compact, inexpensive and retain its heat and mass 
transfer characteristics after extended use. 

Although the contactor may consist of a simple spray with the re- 
sulting air-desiccant interface being the surface of droplets, entrain- 
ment would be excessive since the finely divided liquid is difficult to 
separate from the air. In addition, it would be difficult to remove the 
heat of condensation and dilution. A drop spray onto extended surface 
finned tubes, Figures 7 and 8, prevents atomization of the liquid desic- 
cant and has a better mass transfer capacity per volume than a simple 


spray [33]. 
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Figure 7. Four Row Finned Tube Contactor Dimensions 
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Dimensions for Figure 7 


Toh = 2.00 in 3 = 1,73 in 
Width = 48 in Length = 127 in 
Height = 6.69 in/(four rows) 
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Figure 8. Dimensions for Determining Free Flow Area of Contactor 





32 


The cooling water flows through the tubes in a parallel-counter- 
flow configuration with an even number of tube passes, while the air and 
liquid desiccant flow over the extended surface. Thus, the dehumidifi- 
cation contact surface and cooling coil are conveniently packaged in one 
unit. The cooling water serves to: (1) keep the desiccant cool and thus 
lower the vapor pressure; and (2) cool the air. 

The regenerator contactor has the same physical characteristics 
and performance criteria as the dehumidifier contactor, although the size 
may differ because of different air flow and liquid desiccant recircu- 
lation rates. Therefore, the computer model presented later simulates 
the heat and mass transfer of both the dehumidifier and regenerator con- 
tactors. 

The air flow can be either parallel with or counter to the liquid 
desiccant flow direction, depending on the particular application. Due 
to the dilution of the liquid desiccant as it flows through the contac- 
tor, the exit humidity ratio is not as low as can be achieved with the 
counter flow arrangement. However, the downward parallel flow configu- 
ration does not inhibit spraying the liquid on the contactor and allows 
some separation of the liquid and air by a reversal of flow direction. 

The parallel flow configuration is best suited for the processing of 
large volumes of air, while the counter flow configuration can deliver 
the very dry air needed in manufacturing operations such as the handling 
of hygroscopic materials. The parallel flow configuration represents 


our physical system. 
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Other System Equipment 


Other equipment items, such as the main circulating blowers or 
liquid sorbent pumps, should be selected on the basis of good engineer- 
ing design with cognizance of the difficulty of containing TEG because 


of its propensity to leak past defective fittings or packings. 


LDD System Advantages 


1. Low temperature regeneration as opposed to other thermal cycle air 
conditioners is possible and allows more efficient operation of solar 
collectors. 

2. Regeneration is not very sensitive to heating water temperature, un- 
like other thermal cycle air conditioners. 

3. Solar or recovered waste power can be used to regenerate the liquid 
desiccant. 

4. The system is simple in construction and operation since the cycle 
is without high pressure liquids or gases. 

5. Precise control of air temperature and humidity results from a speci- 
fic TEG concentration and cooling water temperature. 

6. Co-sorption benefits include the removal of bacteria, odors and 
various airborne irritants [1]. 

7. The expense of a separate cooling coil can be saved or a smaller 
after cooler can be specified if chilled water is used in the dehumidi- 
fication contactor. 

8. Excess dehumidification capacity when obtained from low cost waste 
heat or solar energy may be advantageously converted to sensible cooling 


by spray cooling. The spray process cools and humidifies essentially 
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along the lines of constant wet bulb temperature in Figure 4. 

9. In ocean engineering applications dew point temperatures on the 
order of 10 °F below the structure temperature are required to prevent 
condensation [10]. The combination of LDD and cold sea water can produce 
sub-freezing dew point temperatures as the antifreeze qualities of TEG 


prevent contactor frost accumulation. 
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IV. METHOD OF ANALYSIS 


Mathematical Model 

The analysis of the contactor is based on the development of a 
mathematical model suitable for a computer simulation. The substance 
of the model developed is presented here; detailed information may be 
found in the appendices. 

The model of the contactor is a fifth order differential equation 
composed of three first order differential equations for heat transfer 
and two first order differential equations for mass transfer. The 
three heat transfer state variables are the gas, liquid, and cooling 


water temperatures as shown in Figure 9. 
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Figure 9. Locations of Representative Temperatures and Sources 
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TEG Solution thickness 9-007 to 0.010 in 


Fin thickness 0.010 in 
Tube wall thickness 0.035 in 
Tube diameter 0.750 in O.D. 


Dimensions for Figure 9. 


The two mass transfer state variables are the masses of the gas, (water 
vapor plus dry air), and liquid, (water plus TEG). 

All differentials are taken with respect to flow distance in a 
steady state model. This is a reasonable assumption, as start up trans-— 
ients will be negligible within a minute or two of operation from a de- 
termination of the system RC time constants. 

Actual variables used in the program for the mass state variables 
are the superficial mass velocities of the gas and liquid. Thus, the 
five independent variables in this model are ce ty» tes G, and L. All 
the other variables are dependent variables and can be calculated from 
known relationships with knowledge of the five independent variables. 

Two distinct methods were used to develop the mathematical model. 
Appendix A presents a derivation using the traditional method for analyz- 
ing mass transfer operations based on material and energy balances in 
control volumes of a differential section of the contactor. Appendix B 
uses a general method. It develops a bond graph model of the contactor 
from which the governing differential equations rigorously follow. The 


equations developed by the two methods are equivalent. 
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The fiye diffenential equations from Appendix A and B are: 
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Simulation 

These five differential equations and the other supporting equa- 
tions listed in Appendix B which calculate variables such as the mass 
transfer flux, humidity ratio, and TEG weight fraction concentration 
form a subroutine. A master program uses the forth order Runge-Kutta 
method and this subroutine to produce an incremental numerical solution. 

Inputs number 19. The geometry of the contactor is specified with 
eleven inputs. Initial temperatures take three inputs. Flow rates and 


concentrations each take two inputs and a final input controls the in- 
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crement size and the stability of the solution. 

A substantial portion of the subroutine calculates initial values 
of approximately 100 variables including the fluid properties and trans- 
fer parameters discussed in Appendix C. 

A minor portion of the subrountine is repeatedly called during the 
actual solution. Even so, approximately three minutes are required to 
calculate the 2500 increments of a 16 row contactor simulation. Each of 
the five state variables is calculated four times per increment while 
the auxiliary variables are calculated only once per increment. 

The primary outputs of interest are the state variable cS and the 
auxiliary variables: Y, humidity ratio; Po percent TEG weight fraction 
concentration; and Z, the distance in air flow direction. Cross plots 
of s and Y produce operating lines on psychrometric charts while cross 


plots of Po and Z show the regeneration of the TEG. 
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V. RESULTS OF SIMULATION 


Dehumidifier Contactor 

In assessing the performance of the dehumidifier contactor, the 
operating line on a psychrometric chart furnishes the pertinent infor- 
mation: the inlet and exit conditions and the path between them. The 
psychrometric chart, Figure 10a, displays a parametric study of dehumidi- 
fier operating lines for three TEG inlet concentrations, 92%, 94% and 
96% and for three cooling water temperatures, 60°, 70° and 80°F. All 
simulations begin at 92 °F dry bulb and 76 °F wet bulb (approximately 
50% relative humidity) which is a typical outdoor summer temperature 
specification for performance rating of air conditioning units. 

Marks placed on the operating lines of Figure 10a represent equal 
flow distances of four rows (6.69 inches) except for the special mark 
which represents conditions two rows (3.34 inches) from the inlet (top 
of figure). The parameters specified in the figure are the manufacturer's 
recommendations for a 16 row contactor. For all nine simulations the 
mark indicating 8 rows is very close to the 16 row mark (within 2 to 4°F 
and 3 to 4 grains), and represents a 50% savings in contactor acquisition 
cost. Also, a 4 row contactor using 60°F cooling water is nearly equiva- 
lent in performance to a 12 row contactor using 70°F cooling water for 
the same TEG inlet concentration. Such equivalent performance points 
can be the basis for optimizing cooling water operating costs and equip- 
ment acquisition costs for a specific application. 


The dilution vs. flow distance of the TEG solution for the simula- 
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tions of Figure 10a is presented in Figure 10b. In Figure 10b, 

for the same inlet TEG concentrations, the lowest exit concentrations 
correpsond to the lowest temperature cooling water simulations and the 
drier conditioned air. The data on Figures 10a and 10b are from the 
same simulations. 

A correlation of Figures 10a and 10b with the equilibrium humectant 
values of Figure 4 shows that the exit conditions of the 16 row contac- 
tor, operating as recommended by the manufacturer, are at equilibriun. 
For example, from Figure 10b an exit TEG solution concentration of 90% 
corresponds to an inlet concentration of 94% with 60°F cooling water, 

On Figure 10a the 94% simulation with 60°F cooling water has exit con- 
ditions of 60.5°F and 0.005 H.R. or 35 grains. On Figure 4 a dry bulb 
temperature of 60.5°F and a 0.005 H.R. corresponds to a TEG equilibrium 
humectant value of 90% (the Union Carbide data base is used in the compu- 
ter simulation), which is the exit TEG concentration of Figure 10b. 

When psychrometric processes are graphically displayed in the liter- 
ature the path lines are often straight. The path lines for the liquid 
desiccant dehumidification simulation are not straight and can be divided 
into two regions (Refer to Figure 10a.) The first region is nearly all 
latent cooling (isothermal dehumidification) in which the ability of the 
cooling water to perform sensible cooling (dry bulb temperature change) 
is inundated by the exothermic condensation process. In the second reg- 
ion sensible cooling is more pronounced. 

The first region physically represents the flow through approxi- 


mately the first two rows of a 16 row finned tube contactor or approxi- 
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mately 3.34 inches through a contactor 26.75 inches high. The combin- 
ation of a high rate of mass transfer at the inlet (see Figure 11) and 
the associated heat of condensation partially explain the nearly iso- 
thermal region. The exothermic heating due to mass transfer is obvious 
in Figure 12 which shows the variation of the liquid, interface, and air 
temperatures against distance. 

Chemical engineers frequently use a graph of the operating line and 
equilibrium curve in analyzing the absorption of a substance in gas- 
liquid operations. Figure 13 shows these lines for the present system. 
Some simplifications of the calculation procedure can be made if the 
two lines are straight or are straight and parallel [23]. However, for 
the LDD system these lines are not straight because of the temperature 
influence on the equilibrium vapor pressure and hence a computer simu- 
lation is required. 

The operating line of Figure 13 represents the local mole fraction 
of water vapor in the air at the local mole fraction of water in the TEG 
solution. The equilibrium curve represents the equilibrium mole fraction 
water vapor in air at the local mole fraction water in the TEG solution. 

The simulations for Figure 10a were iterated to achieve identical 
inlet and exit TEG solution temperatures. This corresponds to recircu- 
lating the solution through the contactor. Thus, the sump temperature 
for each simulation is slightly less than the air outlet temperature. 
The greatest change in liquid temperature is near the inlet (refer to 
Figure 12) so that the liquid temperature at the 8 row point of a 16 row 


contactor differs little from the liquid temperature at the exit. For 
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this reason the 8 and 12 row points may be accepted as exit conditions 
even thaqugh they are not. Figure 14 verifies this point as the 8 row 
exit marks obtained by iterating the simulation for equal inlet and exit 
liquid temperatures nearly coincide with the 8 row marks of Figure 10a, 
the 16 row contactor. 

Figure 15 shows the results of an air flow rate which is 75% of the 
air flow rate of the Figure 14 simulations. Slightly drier, 1.5 grains, 
and slightly cooler, 1°F, air is the result of a reduced air flow rate 
on the exit conditions. 

Specifying a higher aspect ratio (defined as a greater face area to 
flow distance ratio for a given volume of contactor) will allow the 
treatment of more air. Figures 10a and 14 support this statement. For 
example, the 16 row contactor with a 48" by 127" face area could be 
arranged as an 8 row contactor with a 96" by 127" face area. This high- 
er aspect ratio contactor will process twice as much air assuming con- 
stant TEG recirculation rates and air flow rates per unit face area to 
exit conditions within 2 to 4 °F and within 2 to 4 grains of the 16 row 
contactor exit conditions. 

An indicative result of an 8 row dehumidifier contactor (from Fig- 
ures 4 and 10a) is that with 70°F cooling water and 92% inlet TEG con- 
centration the inlet air at 92°F and 111] grains exits at 72.5°F and 56 


grains or a dew point temperature of 51°F. 


Regenerator Contactor 


In assessing the performance of the regenerator contactor the exit 


TEG concentration is the significant parameter. Figure 16a shows the 
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ability of the regenerator contactor to concentrate 90% and 92% TEC so- 
lution at temperatures of 120°F and 140°F. At constant temperature re- 
generation the exit concentration variations are less than the inlet 
concentration variations. For example, at 140°F regeneration, the in- 
let concentration variation for two simulations is 2% while the exit 
variation is 1/2%. Figure 16b shows the regenerator operating lines on 
a psychrometric chart. Both Figures l6a and 16b are for the same simu- 
lations. 

The regenerator exit conditions are nearly at equilibrium as a 
study of Figures 16a, 16b and 4 will verify. For example, the nearly 
95% exit concentration of Figure 10a for 90% TEG inlet concentration and 
120°F heating water corresponds in Figure 10b to 118°F and 0.021 H.R. or 
147 grains. On Figure 4 a dry bulb temperature of 118°F and 0.021 H.R. 
corresponds to an equilibrium humectant value of about 95% (Union Carbide 
data base). The second equilibrium relationship, equation (2.2), for a 
contact temperature of 118°F and 95% TEG concentration predicts a vapor 
pressure of 24.6 mm Hg which corresponds to 0.0208 H.R. 


24.6 
735.4 





Ze Ce) Uae we = 
Y= 0.622(—-5 z ey 0.622(¢ ) 0.208 (6.1) 


The significant result of the 8 row regenerator contactor simulation 
is that TEG concentrations of 97% to 97.5% can be achieved with 140°F 
heating water, and 95% to 96% TEG concentrations can be achieved with 


120°F heating water. 
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VI. CONCLUSIONS 


1. Contactor equipment costs may be substantially reduced because half- 
sized contactors can provide 95% of the dehumidification and 90% of the 
regenerated TEG concentration as compared to the full sized contactors 
recommended by the equipment manufacturers. 

2. Most effective dehumidification per unit contactor volume occurs at 
the highest ratio of face area to flow distance. 

3. The sensitivity of the regeneration process to the regeneration 
temperature is minimal. 

4. Performance estimates of dehumidification processes which continue 
to equilibrium can be prepared from Figure 4, Triethylene Glycol Equili- 
brium Humectant Values on a Psychrometric Chart. First, assume that the 
exit conditions lie on a vertical line at the cooling water temperature. 
Then, since each point on the line represents a humidity ratio and a 
humectant value, and the sum of the air and liquid masses is constant, 


the equilibrium point can be found. 
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APPENDIX A 


CONTROL VOLUME DERIVATION OF RATE RELATIONSHIPS 


The equations derived in this appendix are concerned with the inter- 
phase transfer of mass (water vapor) and energy which occurs during de- 
humidification of the humid air by the concentrated liquid desiccant. 
Knowledge of the equilibrium characteristics of the system is necessary 
for an understanding of the mass transfer operation and, because of 
simultaneous heat transfer, enthalpy characteristics must also be con- 
sidered. 

Operations such as liquid desiccant dehumidification quickly become 
complex owing to the mass transfer in both the gas and liquid phases and 
the large heat effects which accompany condensation and dilution. Fur- 
thermore, the need to remove the heat of condensation results in a non- 
adiabatic operation. 

Figure 17 represents the direct contact of the gas (dry air and 
water vapor) and liquid (TEG and water) flowing side by side in parallel 
flow and separated by the gas-liquid interface. The changes in tempera- 
ture, concentration, etc., are all differential over the section. 

Based on Figure 17, the balance of this appendix develops: 

(1) An overall material balance for control volume IIT; 

(2) Energy balances for control volumes I, II and III; 

(3) An interface control volume energy balance. 

Next, parameters are separated into catagories and finally, a solution 


sequence is presented. 
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Figure 17. 


Differential Section of Contactor 
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The interfacial surface area in the section is dS. For the finned 
tube contactor the heat transfer area per unit volume, ay is known 
from calculations. Because of the low surface tension of the TEG solu- 
tion complete wetting of the heat transfer area may be assumed so that 
the mass transfer area per unit volume, ay is equal to the heat trans-— 


fer area: 


ay = ay (A.1) 


and, 


dS = adzZ (A. 2) 


where the volume per unit section is dZ. 

Consider the mass transfer operation of Figure 17 as conducted in 
a steady state cocurrent fashion in which Le and G represent the super- 
ficial mass flow rate of the non-diffusing pure TEG and dry air respec- 
tively. Since the values Le and G, are constant through the dehumidifier, 
they are the basis for the material balance. 

At the entrance to the dehumidifier the gas phase contains G total 
lb/ (hr) (sq ft) consisting of non-diffusing Ge, 1b/ (hr) (sq ft) and the 
diffusing water vapor. The water vapor concentration is given by Y the 


mass ratio, or equivalently y the mass fraction. The relationships 


among these variables are: 


(A.3) 


Q 
if p} 
Il 
@ 
o~ 
= 
t 
“< 
w~ 
Il 
a 
+ 
<q 


————— (A.4) 
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y= tty CARS) 


Similarly, the liquid consists of L total 1b/(hr)(sq ft) containing 


x mass fraction of water (or mass ratio X) and L, 1b/(hr)(sq ft) of non- 


S 


diffusing TEG. The relationships among these variables are: 


“ih 

Le = L(1-x) i+x (A.6) 
Popes Lee 

x =>=— (A.7) 
es en 

a ere G2) 


Control Volume III Material Balance 
The quantity of water vapor in the gas is Gy lb/(hr)(sq ft). In 
the liquid it is Lx lb/(hr)(sq ft). Hence the overall material balance 


about control volume III of Figure 17 for the water is: 


iy = Out 
Lx + Gy = (L+dL) (x+dx) + (G+dG) (y+dy) (A.9) 
expanding, 
Lx + Gy = Lx + Ldx + dLx + dLdx + Gy + Gdy + dGy + dGdy (A.10) 


neglecting the second order terms 
d(Lx) = -d(Gy) (A.11) 


or 


Ldx = ~G.d¥ (A.12) 
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The mass transfer rate in lb vapor/(hr)(sq ft) of dehumidifier cross 


section) is [36] 


Nate = ~G.d¥ (A.13) 
where 
— ib water vapor absorbed = ISS ee 
Na (hr) (interfacial surface) earn ley ) (A.14) 


m 


Here y* is the equilibrium mole fraction concentration of water vapor 
over the TEG solution at the interface; oe is the bulk mole fraction 
concentration of water vapor in the gas and Fo is the mass transfer co- 
efficient for large mass transfer rates and is defined by Eq. (C.1). 
Equation (A.14) is a special case for the absorption of one component. 


The gas sensible heat transfer rate across the interface in BTU/ (hr) 


(sq ft) is [36] 


NACA 
eae ~ 1-exp(-N,C,/h,) (Co = oe Gace (ots) 
or 
Iectyt 2 = he(t, - t,)a,dz (A.16) 


where OR = heat capacity of water vapor at constant pressure, 


BTU/ (1b mass) 


Sr 
i 


convection heat transfer coefficient, BTU/(hr) (°F) 

(sq ft interface) and is defined by Eq. (C.2). 

Here he includes the effect of mass transfer on heat transfer by 
accounting for thesensible heat carried by the diffusing water vapor. 


Radiation has been neglected. 
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The sensible heat transfer through the liquid in BTU/(hr)(sq ft) 
is [36] 


Wg Gye > ley Se os et (A.17) 


Energy balances around each of the three envelopes of Figure 17 can 


now be developed. 


Envelope I Energy Balance 


In = Out 


GH. = G. (HH, oF dH.) ae Iegay92 qP NyaydZ(C, (t.-t +A) (A.18) 


where the last term is the enthalpy of the transferred vapor relative 
to a reference state at to where a is the latent heat of vaporization 


at t.. Since 
fo) 


Gey 
Naud = -d(Gy) = -d 1 = ~G,.d¥ (A.19) 
substitution yields 
0 = GodH, + hoay (tg-t, )dZ - Go(C, (tt +A) d¥ (A. 20) 


The enthalpy of the air-water vapor mixture per pound of dry air is 


H, = Cy (tt) + ¥(C,(tp-t +A) (A. 21) 
where Ce = heat capacity of dry air at constant pressure, 
BTU/ (1b) (°F) 
Define 
Cy = C, ae YC, (A.22) 


= humid heat of vapor-air mixture, BTU/(1lb dry air) (°F) 





so that 


He = Co(t, - t) + YA, 


Differentiating equation (A.21) 


dH 


+ = 
G C.dt YC dt, + (C(t, t+ hay 


G A G 
or, 


dH, 


Codt,, ae (C, (tt) + hD dY 
Combining with (A.20) gives the envelope I energy balance 
~G,.C.dt,, = hoa, (to-t,) dz 


Envelope II Energy Balance 


In + heat generation = Out 


4 aes = + 
Ber pa ag §y7t,) + Igy Ant? = Lg CH td) + Q, 


Here 
QQ = Va, (t,-t ,)dZ 
= heat removed by cooling coil, BTU/hr 
where U = over-all heat transfer coefficient between the 


liquid and cooling water, BTU/(hr) (sq ft) (°F) 


cooling heat transfer area per unit volume 


is) 
i] 


sq ft/cu ft 


t. = cooling water temperature, °F 


Substituting equations (A.17) and (A.19) into (A.27) yields 
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(A.23) 


(A.24) 


(A725) 


(A. 26) 


(A.27) 


(A. 28) 
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® - GC, (t,-t d¥ + ha, (t,-t, dz = Lod, + Qe (A.29) 


Here 


Ee) 
o 
W 


N 4,428, = ~GoH dy (A.30) 


the exothermic heat of dilution (mixing), BTU/hr 
where Hy = heat of dilution, BTU/(1b water transferred) 
Substituting equations (A.28) and (A.30) into (A.29) yields 


~G, (H)tC,, (t,-t) d¥ + hy ay, (t,-t, dz = Lo dH, Ua, (t,-t.)dzZ (Aza) 


The enthalpy of the TEG - water solution per pound of TEG is 


Hy C.(t, -t) + XC, Ct, -t) {A.32) 


where C. = heat capacity of TEG at constant pressure 


BTU/ (1b) (°F) 


Car = heat capacity of liquid water at constant pressure 


BTU/ (1b) (°F) 
Define 


Coy = C. + XC. (A.33) 


= the humid heat of the TEG-water solution, 


BTU/ (1b TEG) (°F) 


Differentiating equation (A.32) 


= of = = 
dH (C, XC, dt, + one (t, t dx (A. 34) 





or, dH, = Co pot, = Car (t)-t dx 
Using equation (A.12) 
= G —_ 
dH, = Cop dt, - = aL (t, t dy 
S 


Combining with (A.31) gives the envelope II energy balance 


et. dt = (h,a,d2Z-G.C apd) (ts -t 1 7Ss H_dY-Ua,(t,- 


5 SL L 8 D C°L ee 


Envelope III Energy Balance 


The overall enthalpy balance for control volume III is 


In Out 


LH + CoH, + ors = L, (A, + dH, ) + Go (CH, + dH.) + Q. 


Rearranging 


he > Ee ees 


Substituting into equation (A.37) 


ec. dt -G gle. dt, + (C, (t, -t Sle Cc 


at LC (t 


~t +A tH) d¥] - 


AL’ L 


UVa, (t, - t.)dz 


Interface Control Volume Energy Balance 


In Out 


Ny aydZQ +C, (tt )) + Iggy 42 = Na,92(C, 7 (t-t)) + Io 292 
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(A.35) 


(A. 36) 


(Anen) 


(A. 38) 


(A.39) 


(A.40) 


(A.41) 





Rearranging, 
Mee Co AL into) * %o? = MHP, fFi 4) -hg(t,-¢,)) 
Solving for t. 
«= Nay CC, (tt tC, t tA) ta, (ht, thot.) 
SeMaug? Hou c: 
Parameter Classification 


1. Known quantities for which inlet values are given 


to 3eR Ho enue ys 
the x, Hy > bx, 
a 

2. Constants 
Go» Cas Che aw? 
Los Car? Cos Hp» 
1005 to 


3. Variables with known functions 


* 
ayy aye sat ol fo hee hy > is Me 
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(A. 42) 


(A. 43) 


By starting at the top of the dehumidifier with the known quantities 


and solving the equations at each incremental step, dZ, 
quantities will vary; 


crag eHencany, Co 


peeceie bx) C_, 


the following 
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* 
Ci Ny» Qn? y ’ Dee 
Eee Qe 
and will be obtained after each step and at the exit of the dehumidifier. 


Solution Sequence 
A solution sequence based on the differential section of Figure 
17 is: 


aoe 
Vege (A.44) 


water mass ratio, 


x 





xX = jee (A.45) 

humid heat, 

Cy = Cc, oF YC, (A. 46) 
liquid humid heat, 

Cor = C, + XCar (A.47) 
gas enthalpy, 

Ho = Co(ta-t)) + YA, (A.48) 
liquid enthalpy, 

HL 3° Co, (t,t) (A.49) 


interface temperature, t.> from Eq. (A.43);3 
equilibrium vapor pressure above liquid, p*, from Eq. (2.2); 
mole fraction concentration in the gas in equilibrium with the bulk 


liquid from Henry's Law, 


se pe (A.50) 





mole fraction concentration in gas, 


a) 
Ym 0.622 + 0.378 y 


mass transfer flux, 


= eyes 
Ny Fo aie ) 
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(ASS) 


(A.52) 


convection gas heat transfer coefficient corrected for simultaneous 


mass transfer 


De Naa 
1-exp(-N,C,/h)) 


(A.53) 


This completes the supporting equations for the solution sequence. The 


derivatives for the solution follow: 


t 
dt, = ~h,a, (tg-t,) 


dz eae 
ie es = Z 
= Tice Ne + Cy, Gy gz) (t,-t,) +H, (-6 


Va, (t,-t,)) 





dt. = Ua_(t.-t.) 
aia Oe c Le 

f “AL 
dy oN 
dz “4 ay 

S 
oe: dy. 
az 1 Gee dz 

S 

dy_ 1 wy 
dz ae. y)2 dZ 
cpr Lace 
dz dz 


(1 + x)? 


a¥y_ 
SadZ 


(A.54) 


(A.55) 


(A.56) 


(A.>7) 


(A.58) 


(A.59) 


(.60) 
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ay) = G, oe (A.61) 
atx) = -c, & (A.62) 


Before simulation of the contactor with the solution sequence can 
proceed, (1) input values must be initialized to correspond to the phy- 
sical system; (2) equations for fluid properties as functions of temper- 
ature and concentration are needed; (3) heat and mass transfer coeffici- 
ents for the physical system must be predicted or measured; (4) a ra- 
tional selection of the minimum number of variables necessary to specify 
the state of the model must be made. 

Equations (A.54) through (A.64) above are not all independent. 
Appendix B describes through bond graph techniques the selection of the 
minimum number of state variables to represent the LDD contactor. The 
results of Appendix B indicate that differential equationa (A.54), (A.55), 
(A.56), (A.63), (A.64) above, plus necessary auxiliary equations, will 


model the contactor. 
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APPENDIX B 
BOND GRAPH CONTACTOR MODEL 
Introduction 

The bond graph method developed by Professor H.M. Paynter of MIT 
for systematic derivation of differential equations for the transport of 
power is described briefly: Bond graphs are a powerful analytical tool 
as they give a physical presentation of power flows by maintaining the 
topological structure of the system while allowing separate determin- 
ation of the sign of positive power flow and distinction between depen- 
dent and independent variables. While each engineering discipline has 
its own calculation methods, bond graphing is compatible with all dis- 
ciplines and is invaluable for model development at interfaces, i.e. 
electrical and fluid. Finally, the transformation from a bond graph to 
equivalent differential equations follows an exact method. 

Bond graph theory has been treated rather thoroughly elsewhere [13], 
[14], [25], [34], [35]. Thus, only those aspects necessary for under- 
standing the present model development will be discussed and only a 
meager sampling of the analytical potential of bond graph theory is pre- 


sented. 


General Bond Graph Model 

A bond graph representation of a control volume around the gas- 
liquid interface is shown in Figure 18. This bond graph is a distributed 
resistance-capacitance model of the interface and can be replicated n 
times to provide a model of any desired accuracy. Dynamic simulation in 


time and space is possible with this general model. 
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Figure 18. Bond Graph Representation of Coupled Thermal and Mass Transfer 
Operations in Liquid Desiccant Dehumidification 
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The zeros represent parallel junctions at which the sum of the 
heat or mass flows equal zero. A distinct temperature or chemical po- 
tential is associated with each zero junction. Two types of capaci- 
tances are represented on the graph: thermal, C. and mass, ae There 
are two types of thermal resistance: convection resistance in the 
fluids, Re and interfacial resistance, R,- Mass transfer resistance 
is also of two varieties: diffusion resistance, Ry and reaction re- 
sistance, Ro 

The dotted line is a reaction coupling from the mass transfer por- 
tion of the bond graph to the heat transfer portion of the bond graph 


and represents the heat of condensation and dilution. 


Specific Bond Graph Model 


Looking at a time invarient solution of the general model allows 
the n replications to be done with a dynamic simulation routine if the 
differentials are taken with respect to dZ. Figure 19 shows this speci- 
fic pineiditied model with the addition of a cooling medium with heat 


capacity C,. and temperature t 


ft £ 

Bond graphing allows a determination of the number of independent 
variables by counting the number of storage elements in the model. The 
temperatures oe t, and te are chosen to represent the mass transfer 
independent variables because of their physical significance and their 
association with the storage elements. Figure 20 is a word representation 


of Figure 19 and the contactor. (Note that the t's represent temperature 


differences.) 





73 


suotjeisdQg Avjsueal sseW pue [eWAsYyL JUeOTJTUSTS Jo ydeay puog AoJoePJUOD “GT aAN3TY 


'y 








cm~e 
=e 











[=] oe ee 


ray h 


co —_o— 
= 
Bas (=) 
= 
Ae 








74 





At At At ‘ 
air interface ———— desiccant ——2— cooling 
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Figure 20. Word Bond Graph of Contactor 


Bond Graph Symbols 


On Figure 19 the bonds with arrows represent signal flows with no 
feedback along the same bond. For example, the interfacial resistance 
to mass transfer, R.> is dependent on the interface temperature. Solid 
bonds represent the power flow associated with heat or mass transfer. 
Ones represent series junctions at which the sum of the efforts, either 
temperature or chemical potential, is zero. 

In a true bond graph the product of the effort variable (above the 
bond) and the flow variable (below the bond) is power. Figure 19 is a 
pseudo bond graph because the concentrations are used as the mass trans- 
fer efforts rather than chemical potential and because heat flux is used 
as the thermal flow rather than entropy. With a pseudo bond graph the 
reaction couplings must be scrutinized as standard bond graph rules do 
not hold, but with consistent units, pseudo bond graphs follow the re- 
mainder of the rules for regular bond graphs. The pseudo bond graph was 


selected for the physical significance of the effort and flow variables.[14] 
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A "T" in a signal bond indicates the signal, t, or N,, is trans- 
formed before use. The dotted line indicates the transformer modulus 
is controlled by oe Linearity is not implied by the bond graph. Two 
mass flux controlled heat sources, Sey and Seq represent the latent heat 
of condensation and the heat of dilution. 

Heat flux in BTU/(hr)(sq ft) and mass flux in lb water/(hr) (sq ft) 


are the flow variables used in Figure 19. Table 1 lists the units of 


the thermal bond graph. 


Variable Units 
Effort °F 
Flow BTU 
(hr) (sq ft) 
Resistance (hr) (sq ft) (°F) 
BIU 
Capacitance BTU 


(hr) (ft) (°F) 
Table 1. Thermal Bond Graph Units 


Dimensionally, the product of resistance and capacitance is (ft) or 


the same as the solution increment. 


Differential Equation Model 


The differential equations which result from a rigorous application 
of bond graph techniques have been presented as equations (4.1) through 
(4.5) and are equivalent to equations derived in Appendix A. The com- 
puter simulation uses: (1) differential equations (4.1) through (4.5); 


(2) the auxiliary equation for tis Eq. (4.6); (3) the equilibrium relation- 
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ship for vapor pressure over TEG, Eq. (2.2); and (4) the following se- 


quence of auxiliary equations: 


mass fraction in liquid, 


be 
= 


L 


percent weight fraction concentration of TEG, 
Po = 100 (1 - x) 

interface temperature from Eq. (4.6), 
ti. Bee NA 

equilibrium vapor pressure over TEG from Eq. 


k= 
P fet(t, »P,) 


(B.1) 


(B. 2) 


Henry's Law for the mole fraction concentration in the gas in equilibrium 


with the bulk liquid, 


p* 


oe 


y* = 


mass fraction in gas, 


mole fraction concentration in gas, 


ie y 
yee 
0.622 + 0.378y 


mass transfer flux (Fo is derived in Appendix C), 


(B.3) 


(B.4) 


(B.5) 


(B.6) 
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heat of dilution taken as 10% of the heat of condensation and sensible 


heat associated with the transfered mass, 
Q%® = Ny (100 + (t, - t,)) (B.7) 
mole fraction concentration in liquid, 


——— (B.8) 
“a gua9h = 7.334 x 


and humidity ratio, 


ss (B.9) 


The computer model is now complete. 
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APPENDIX C 


FLUID PROPERTIES AND TRANSFER PARAMETERS 


This appendix presents the essential but less than obvious relation- 
ships. The development of an equation for the vapor pressure over TEG 
from graphical data is discussed in Chapter II and is representative of 
the method of the numerous curve fits presented here. The remaining de- 
tails may be found in the computer program listing. Appendix E. 

Essential to the model is a numerical value for the mass transfer co- 


efficient, F This coefficient is predicted from the more easily ob- 


G? 
tained gas convection heat transfer coefficient by means of the well 
accepted Chilton-Colburn analogy. Agreement of the analogy with newer 


data on mass transfer is very good [32]. The mass and gas heat transfer 


coefficients are related by 





h 0.667 
ae elas 
lee eee Cy 
B 
Fo = mass transfer coefficient, (1b water)/(hr) (sq ft) 
where CB = heat capacity of air at constant pressure, 


BTU/ (1b mass) 
Pr = Prandtl number 
Sc = Schmidt number 
h, = heat transfer coefficient between the gas and 


laminar TEG film, BTU/(hr) (sq ft) (°F) 


A numerical value for he is based on an empirical correlation for a 
representative finned tube extended surface heat exchanger reported by 


Rohsenow [30]. 
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0.681 O335) 2072 270.113 


x k 
ho = Oa B (tomax ) We) @) (2) (ee) 
D u . 
ie B 
where 
D = root diameter of tube, ft 


k_ = thermal conductivity of air, BTU/(hr) (ft) (°F) 


G = mass flow rate at minimum cross section, 
1lb/ (hr) (sq ft) 
s = distance between adjacent fins, in 
1 = fin height or half the mean fin distance 
between tubes, in 
t = fin thickness, in 


UH, = viscosity of air at bulk temperature, lb/ (hr) (ft) 


Eq. (C.2) shows excellent agreement with Kays [15] plate finned circu- 
lar tube data but which geometry differs from the modeled contactor. 


A curve fit developed for the Prandtl number for air is 
Pr = 0.721 exp(-196E-6 t,) (C.3) 
To compute the Schmidt number 


U 
B 
(C.4) 
p DL 





Sc = 


values for the viscosity, density and diffusivity are needed. 


Curve fits led to 


0.159 
Ug = 0.0223¢t,) (c.5) 


359 (460 + t.) 
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while the diffusivity is [2] 


TONS Ge Cees 


ep ey oe ee 
DY e+ 901 (C.7) 


The thickness of a falling film of the water TEG solution on 


the vertical aluminum fin is [26] 





3 ul 0.333 
6 7) (c.8) 
prs 
where 6 = film thickness, ft 


Up = film liquid viscosity, 1lb/(hr) (ft) 


. = liquid loading per unit width of fin, 1lb/(hr) (ft) 
Py = density of film liquid, lb/(cu ft) 
g = acceleration due to gravity, eC/ne 


The mass flow rate per unit width [ can be calculated from the liquid 
recirculation rate but values for Up and PD for TEG-water solutions 


are obtained from curve fits 


(1.487 - 0.727 In P_) 


By? ty iT (c.9) 


4.46 E-7(P 


Up 


(C. 10) 


54.88 + 3.73 In Po - 0.0262 ty 


Pp 


The liquid phase heat transfer coefficient was approximated as 


conduction dominated between the air-liquid interface and the liquid 


mid-point by 
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= 2 (Car) 


from the curve fit for the TEG-water solution conductivity 
ky = 0.789 - 0.14 In Pr - 131E-6 ty (C.12) 


The heat transfer coefficient from the liquid mid-point to the 
aluminum fin was approximated by the empirical correlation for a 


falling film [17], [26] 


0.333 
ate 0.111 
Des 4. 
i= 7? 0.333 7 (C.13) 
il! 10s ena D 
D 
for 
= 7100 
BD 


where L = length of heat transfer surface, ft 


°D 


Cor defined in Eq. (A.33) 
A curve fit for the specific heat of the TEG-water solution is 


Cy = 0.4977 + 0.000444 ty oF XC ar, (c.14) 


A cautionary note with Eq. (C.13), [26] suggested a range of validity 

on L as 0.4 to 6 ft. Accordingly, the contactor height 1.12 to 2.24 

ft was used for L. This gave values for hy from 40 to 60 while values 

for hy ranged from 150 to 200. These values indicate that Eq. (C.13) 

may not be used since the heat transfer coefficient can not be less 

than that of pure conduction. Exclusive use of the conduction coefficient 


hy is correct. The best approach is to calculate both h, and hy and use 


W 
the larger for the entire liquid film. Sole use of hy produces slightly 
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cooler operating lines as shown by Figure 21. In addition the typical 
overall heat transfer coefficient from the liquid to the cooling water 
increases from 14 to 20 while the typical overall heat transfer 
coefficient from the air to the cooling water increases from 6 to 7. 
The comparison in Figure 21 shows that the simulation results are not 
materially affected. 

A curve fit was developed from tube side water heat transfer co- 
efficient data [17] 


Do = 26 0.8 


= = gules enh 
h, = (0.905 - 0.1863 In ( 0.0833 )) (163 +1.731 t,) Ve (c.15) 


where 
o. = tube wall thickness, ft 


Ve = water velocity through tube, ft/sec 


The free flow area of Figure 8 is slightly reduced by the TEG-water 


solution film thickness so that the ratio of free flow area to face 


area is 
os a - 
Gc au NG SA a cae (C.16) 
T 
sh 


where P = number of fins per inch, 1l/in 


poh 


horizontal tube spacing, in 


The overall heat transfer coefficient from the gas to the cooling 


water is calculated by 








il 

v= 1 a5, 8y (C.17) 
i jhe eaie 4. 
mG L t w t ig 
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where 
2 3 
a. = average tube surface area per unit volume, ft /ft 
k = conductivity of tube wall, BTU/ (hr) (ft) (°F) 
5 = thickness of tube wall, ft 


The calculations were done with the fin efficiency applied to both 
the fin area and the tube area. 
Applying the fin efficiency to the total area instead of just the 


fin area accounts for 97% of the heat transfer: 


. 66 


TECDarsOy 90 os (C.18) 


i: 
ly ps 
“Hy 


af [4 | 


The remainder of the contactor geometry and transfer parameter 
calculations is straight forward. A complete listing may be found in 


Appendix E. 
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APPENDIX D 


TRIETHYLENE GLYCOL PROPERTIES 


Triethylene glycol is a colorless, combustible, hygroscopic water- 
soluble liquid used as a chemical intermediate, solvent, bactericide, 


humectant, and fungicide. 


Physical Properties [12], [28]: 


Formula CH, -0- CH, - CH, - OH 


| 
CH, - 0 - CH, ~ CH, - OH 


Molecular Weight 150.2 
Boiling Point at 1 atm 545.9 °F 
Vapor Pressure at 77°F <0.01 mm Hg 


Specific Gravity at 77°F 1.119 g/ee 


Freezing point 19 °F 

Pour Point -73 °F 
Viscosity at 77°F 37.3 centipoises 
Surface Tension at 77°F 45 dynes/cm 
Specific Heat at 77°F 0.53 BTU/ (1b) (°F) 
Flash Point 320: °F 

Fire Point 330 °F 


Conductivity at 77°F 0.13 BTU/ (hr) (ft) (°F) 
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Desired Desiccant Properties [33] 


A liquid desiccant selected for a dehumidification system should 


have the following characteristics: 


Le. 


Suitable vapor pressure characteristics; should be capable of con- 
centration at low temperatures 

Noncorrosive 

Odorless 

Chemical stability; should not break down or vaporize in the range 
of interest 

Low viscosity and good heat transfer characteristics; must not 
crystallize or solidify near the operating range 

Nontoxic and nonflammable 


Widely available and inexpensive 


Triethylene Glycol Desiccant Properties [12], [28] 


ile 


TEG has a suitable vapor pressure curve and can be concentrated at 
temperatures as low as 120 °F 

TEG is compatible with steel if slight TEG discoloration due to iron 
contamination is acceptable. Otherwise stainless steel, aluminun, 

or resin coatings should be used. 

TEG is practically odorless. 

TEG is a superior choice considering both its hygroscopicity (item 1) 
and its high boiling point. 

TEG has a low viscosity and good heat transfer characteristics. 


TEG is low in chronic oral toxicity. It has been fed in the drinking 
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water of rats over a two year period at a 4% concentration with no 
apparent ill effect. Prolonged inhalation of saturated vapors of 
TEG has produced no ill effects on humans. Flash and fire points 
determined by the ASTM Cleveland Open Cup method for pure TEG are 
320 °F and 330 °F respectively. TEG presents a minimal fire hazard 
in storage or handling. 

TEG is a byproduct of automotive antifreeze. The cost of TEG com- 


pares to the cost of antifreeze. 





APPENDIX E 


COMPUTER PROGRAM LISTING 


AND 


SAMPLE OUTPUT 
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